OPTIMAL CONTROL OF A FREE BOUNDARY PROBLEM: 
ANALYSIS WITH SECOND ORDER SUFFICIENT CONDITIONS 



HARBIR ANTIL*, RICARDO H. NOCHETTQt, AND PATRICK SODRE* 

Abstract. We consider a PDE-constrained optimization problem governed by a free boundary 
problem. The state system is based on coupling the Laplace equation in the bulk with a Young- 
Laplace equation on the free boundary to account for surface tension, as proposed by P. Saavedra 
and L. R. Scott 16 . This amounts to solving a second order system both in the bulk and on the 
interface. Our analysis hinges on a convex constraint on the control such that the state constraints 
are always satisfied. Using only first order regularity we show that the control to state operator is 
twice Frechet differentiable. We improve slightly the regularity of the state variables and exploit this 
to show existence of a control together with second order sufficient optimality conditions. 
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1. Introduction. Free boundary problems (FBPs) are challenging due to their 
highly nonlinear nature. Besides the state variables, the domain is also an unknown. 
FBP find a wide range of applications from phase separation (Stefan problem, Cahn- 
Hilliard), shape optimization (minimal surface area), optimal control problems with 
state constraints, fluid dynamics (flow in porous media), crystal growth, biomem- 
branes, electrowetting on dielectric, to finance. For many of these problems there is a 
close interplay between the surface tension and the curvature of the interface pTl \n\ . 
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Fig. 1.1. fl-y denotes a physical domain with boundary dQ-y = S U F-y. Here S includes the 
lateral and the bottom boundary and is assumed to be fixed. Furthermore, the top boundary F-y ( dotted 
line) is "free" and is assumed to be a graph of the form {xi,l + ^{xi)), where 'y £ (0,1) denotes 
a parametrization. F-y is further mapped to a fixed boundary F = (0, 1) and in turn the physical 
domain fi-y is mapped to a reference domain Q, = (0, 1)^, where all computations are carried out. 

Of particular interest to us is the control of a model FBP previously studied by P. 
Saavedra and L. R. Scott in [T2] and formulated in graph form; see Figure [T] where 
the free boundary r..y is the dotted line. The state equations 
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1.2b I involve a Laplace 



equation in the bulk and a Young-Laplace equation on the free boundary to account 
for surface tension. This amounts to solving a second order system both in the bulk 
and on the interface. Below we give a detailed description of the problem. 

Let 7 € (0; 1) denote a parameterization of the top boundary (see Figure [l]) 
of the physical domain fl^ C ft* C with boundary dil^ := U S, defined as 

f}* = (0,l)x(0,2), 

= {{xi,X2) : < xi < 1, < a;2 < 1+7(2:1)} , 
= {{xi,X2) : < xi < 1, X2 = 1 +7(a;i)} , 

r = {{xi,X2) : < xi < 1, 2:2 = 1} . 

Here, fl* and E are fixed while flj and deform according to 7. 

We want to find an optimal control u € Uad C (F) so that the solution pair 
(7, y) of the FBP approximates a given boundary 7^ : F — )• M and potential yd'-^* ^ 
]R. This amounts to solving the minimization problem 



1 2 1 

ToinJ {j,y,u) := 2 IIt - 7<i|lL2(r) + 2^\y ' 

subject to the state equations 
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(1.2b) 



the state constraints 



(1.2c) 



|da:i7(xi)| < 1 a.e. xi € F, 

with dx^ being the total derivative with respect to xi, and the control constraint 

u £ Uad (1.2d) 

dictated by Uad, a closed ball in (F). Here > is the stabilization parameter; v 
is a given data which in principle could act as a Dirichlet boundary control; 



n [7] d. 



is the curvature of 7; and k > plays the role of surface tension coefficient. 

Optimal control of partial differential equations (PDEs) allows us to achieve a 
specific goal ( 1.2a) with PDE ( 1.2b I and other constraints ( 1.2c[ )-( 1.2d I being satisfied 
and can be highly beneficial in practice (see [19] for more details) . For example using 
the reverse electrowetting i.e., by applying a control to change the shape of fluid 
droplets, one can generate enough power to charge a cellphone |13| . by mere stroll in 
the park. There has been various attempts to solve optimal control problems with 
a FBP constraint. We refer to (TUl E] for control of a two phase Stefan problem in 
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graph formulation and 4 for the same problem in level set formulation. Paper [14] 



discusses optimal control of a FBP with Stokes flow. Even though problem (1.2a) 



(1.2d) is relatively simple, it captures the essential features associated with surface 
tension effects found in more complex systems, and allows us to develop a complete 
second order analysis, based on [19], which is absent in the existing literature on FBP. 

There are several methodologies to formulate FBP depending on the role of the 
free boundary. We deal with a sharp interface method, written in graph form (see 
Figure [l]), for which the interface is governed by the explicit elliptic PDE 

-kV. [7] + duV ^ u. 

A similar approach was used in |10[ lllj for the optimal control of a Stefan problem, 
but without the full accompanying theory developed here. Alternative approaches to 
treat FBP are the level set method and the diffuse interface method [HI H] ■ 

We use a fixed domain approach to solve the optimal control free boundary prob- 
lem (OC-FBP). In fact, we transform fl^ to Q ^ (0, 1)^ and r..^ to F = (0, 1) x {1} 
(see Figure [l]), at the expense of having a governing PDE with rough coefhcients. 
This avoids dealing with shape sensitivity analysis [T71 [?]• We refer to for a 
comparison between these approaches applied to a FBP. Using operator interpolation 
[18j we demonstrate how to improve the existing regularity of state variables derived 
earlier in jl6j . which turns out to be instrumental to derive the second order sufficient 
condition. One of the challenges of an OC-FBP is dealing with possible topological 
changes of the domain by introducing state constraints. Our analysis provides con- 
trol constraints which always enforce the state constraints i.e., we can simply treat 
OC-FBP as a control constrained problem. 

We have organized this paper as follows. A detailed problem description on a fixed 
domain is given in [section "2] We introduce the Lagrangian functional to formally 
derive the first order necessary optimality conditions in [section "3] We present a 
rigorous justification of the Lagrangian results in the remaining sections. To this end. 



we introduce a control to state operator in section 4 and show that for a particular set 



of admissible controls it is twice Frechet differentiable. Finally, we write the optimal 
control problem in its reduced form and show the existence of a control under slightly 
higher regularity together with second order sufficient conditions in [section 5[ 

2. OC-FBP on Reference Domain. For simplicity we consider the FBP 



(1.2b) with linearized curvature. To analyze the minimization problem (1.2), we 



map the physical domain 57^ onto the fixed reference domain il = (0, 1)^. This results 
in an optimal control problem subject to PDE constraints with nonlinear coefficients 
depending on 7 but without an explicit interface. The idea is to map the unknown 
domain onto the fixed domain fl = (0, 1)^ using the inverse of the Lipschitz map, 
: — >• ri^, defined as 

*(a;i,X2) = (xi, (1 -t-7(xi)) e f^T,, for (xi, X2) G f^. (2.1) 



Since 7 is Lipschitz continuous according to the state constraint (1.2c), we deduce 
that I7I < 1/2 which in turn implies that 5* is invertible. Furthermore, the inverse of 
is also Lipschitz. Moreover, it becomes routine to check that the Laplace equation 
Ay = in fi^ and d^y on can be written as 



div (A [7] Vj/) = inn, A^Vy •z/(l+|d,,7|') 



-1/2 

on F, 



where v = [0, 1]^, and A : (T) L°°{n) is the Nemytskii operator \l9i Chapter 
4] defined by 



^[7] = 



1 + 7(2:1) 



cia;i7ia;ija;2 1+7(2:1) 



(2.2) 



-1/2 



It is convenient to write A2 2 [7] as $ [7] := f {j{xi),dxij{xi)x2) , where ip{a,b) 
(l + 52)/(l + a). 

To simphfy the exposition we make the foUowing assumptions: 

(Ai) Linearized curvature: Hun [7] — d^^7 I 1 +|d2:i7| 

(A2) Scaled control: u ~ u (^1 + \dj:^^f^ 

(A3) Dirichlet boundary condition: we identify v on the boundary with the trace 
of a function v G Wp {fl). 
These assumptions are not crucial. We refer to jTHi Section 4] for the FBP involving H 
instead of Hun ■ The scaling of the control avoids unnecessarily complicating the right 

hand side of(2.3b| below, which would contain u 1 1 +|da;i7| ) instead of simply u. 
Finally, the lifting of the Dirichlet v data is standard practice. 

Under these assumptions and the application of the map the optimal control 
problem (1.2) becomes 



■ ^ / ^ 1 II II 2 I '^a II , II 2 , II ||2 

subject to the state equations (7, y) e VF^ (F) x Wp (fl) 

( - div {A [7] V (y + w)) = mil 

1 — '«d^^7 + A [7] V (y + w) • = u in F 



(2.3a) 



(2.3b) 



together with the state constraint (1.2c I and the control constraint (1.2d). 

We point out that the part of the cost functional J' in (2.3a I corresponding to 
state y is defined over the fixed domain instead of il^ without taking into account 
the determinant of the Jacobian iD^t'l of because 



yd\ dx^ 



ly + v-ydHB^S/ldx 



-7|dx, 



whence 



P 2 /" 2 ^ /" 2 

\u + V ~ yd\ dx < \y + v-yd\ dx^ < - / \y + v - yd\ dx, 
Jn Jo.-, ^ Jn 



where we have used I7I < 5. Therefore, we define over without |D5'|. 

In order to derive the first and second order optimality conditions in later sections, 
we need to compute the first and second order directional derivatives of A, which 
in turn requires computing the directional derivative of the Nemytskii operator ^ 
defined above. To simplify notation, we drop the evaluation of 7 and dx-^j at xi. The 
derivative of $ in the direction h at (7, X2 dx^^) is given by 



D$[7] (h) = daV{l,X2 dx^j)h + dbip{'y,X2 dx^j) X2dx^h 
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where 



l+b'' 2b 
"(l+a)2 1+a 



Furthermore, we obtain the following representation for DA in terms of h and h' 
DA [7] (h) ■.= Ai h]h + A2 h]d,,h 



1 



-X2 

-X2 X2dbip{'J,X2dx^j) 



whence the remainder ^Ra [7, h] at 7 in the direction h reads 

D\a [7, h]:^A[j + h]~A [7] - DA [7] (h) 

and 

h,h]\ 



d.,h, (2.4) 



(2.5a) 



The Hessian of (p is 

VV(a,fe) 



hm 



dl'p(a,b) dabV{a,b) 
dbaf{a,b) dlip{a,b) 



L~(r) 



= 0. 



(2.5b) 



-fc 

(l+a)3 (l+a)2 

-b 1 

(l+a)2 1+a 



The second order derivative of $ in the direction hi followed by /12 evaluated at 
ij,X2d^,-r) is 

, X2 drc-^j) /l2^1 + (?a6<y3 (7, X2 d^j^j) X2h2 d^^hi 

+ dab^P (7: ^2 dxij) X2 dx^h2hi + dlifi (7, X2 d^^l) X2 dxi h2 djjj hi . 
Finally, we obtain the following representation for if A in terms of hi and /i2 



I^A[7] {h2,hi) = 





D^$[7](/i2,/ii) 



whence the remainder IHda [7, ^i, ^2] at 7 reads 

5Hda [7, hi,h2] ■■= BA [7 + h2] (hi) - BA [7] (hi) - T^A [7] (/ii, /is), 

and 

||9iDA [7, /ll,/i2] 



lim 



iL°°(a) 



l''ilw^(r)-+0 l^llH'4(r)l^2|vi/^(r) 



0. 



(2.6) 



(2.7a) 



(2.7b) 



Proposition 2.1 (bounds on A). Applying the state constraint ( |l.2c[ ) to ( |2.2[ ), 

positive ct 

DA [7] (/i)| 



(2.4) and (2.61 implies the existence of a positive constant Ca < 00 such that 



lL~(n) 



sup 
sup 



lL°°(n) 



D?A[7] (/ii,/i2; 



L°°(n) 



(2.8) 
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3. Lagrangian Formulation. In this section we formally derive, ignoring the 
state constraint (1.2c I, the first order necessary optimality conditions using the La- 
grangian approach described in [19] . For a rigorous analysis of the existence of La- 
grange multipliers in Banach spaces we refer to |23j . It is well known that for a convex 
optimal control problem with linear constraints, the first order necessary optimality 
conditions are also sufficient conditions [19, Lemma 2.21]. However, despite lineariz- 
ing the curvature via assumption (Ai), the state equations (2.3bl are still highly 
nonlinear, and the optimization nonconvex. We will derive the second order sufficient 



optimality conditions in [section 5 



Let s, r denote the adjoint variables corresponding to states 7, y respectively, and 
p, q be conjugate indices i.e., l/p+l/q = 1 withp > 2. Then the Lagrangian functional 
is given by 



^l,y,u,r,s) -.^ J {'^,y,u) + / div (A [7] V (y + u)) rdx 

Jn 

+ / (^Kd^_^"f — A['y]V {y + v) ■ 1/ + uj sda. 



(3.1) 



Additionally, if {'y,y,u,r, s) is a critical point for the Lagrangian £, then the first 
order necessary optimality conditions are 

(B{^,y,r,s}C{j,y,u,r,s),h)=0 yh e [w^, {T),W^ {n),W^ {n),Wl{T)} , (3.2a) 
(Di,£(7, y,u,f, s), u - w) > yueUad, (3.2b) 

where {} denotes a list, e.g. (D^ £(7, y, u, f, s), /i) = for all h e (F). Therefore, 
computing {'j,y,u,f,s) requires solving the nonlinear system (3.2). In practice this 
can be realized using techniques described in [3l [121 IE] ■ To solve variational inequal- 
ities of first and second kind we refer to [3] for relaxation and augmented Lagrangian 
techniques and to [S] for semi-smooth Newton methods, and references therein. The 
remainder of this section is devoted to the derivation of the equations satisfied by 
(7, y, u, f, s) using the nonlinear system above. 

Since ^D^-j £(7, y, m, f, s), ft.) — implies that (7,?/) solves the state equations 



(2.3b), we focus on the adjoint equations {^{^^y} £(7, 2/, u, r, s),hj — 0. Using Green's 
theorem and assuming smoothness, the Lagrangian C can be rewritten as: 



^y,l,u,r,s) ^ J {^,y,u) + / div (A [7] Vr) (y + dx 

Jn 

+ / {rA [7] V {y + v) - {y + v) A [7] Vr) • ly dcr 
Jdn 

+ / (^K"fdl^s-sA['y]S/{y + v)-iy + us^da + K{dxj^"fs-jdx-,s) 



(3.3) 



Imposing (Dj, £(7, y, u, r, s), /i) = to (3.31 implies that for every h e Wp (il) 
div {A [7] Vf) hdx — Xy (y + v ~ yd) h dx 



+ rA['y]'Vh ■ lyda - sA [7] Wh ■ i^da 



an 



(3.4a) 



Next, without loss of generality {C^ (ft) is dense in Wp (fi)), we obtain 

div {A [7] Vf)hdx = Xy [ (y + v- ya) hdx \fh € {Q) , (3.4b) 



whereas, using that A [7] V/i • v can be chosen arbitrarily on dfl we deduce from ( 3.4 1) 
and K^) that 



f-s|r = 0, f|s-0. 



(3.4c) 



In view of (3.4d-c), the strong form of the boundary value problem for r is: seek 
f e (n) such that 



- div (A [7] r) = Xy{y + v - yd) in Q 
f — s on r 

f = on S. 



(3.5) 



Next we employ the same technique to obtain the equations for the second adjoint 
variable s: we impose (iXf £(7, y, u, f, s), ft.) =0 to (3.3 1 and make use of the boundary 
conditions in (3.5) to obtain for every h G (F) 

- / Kdl^shda= / (7-7d)/idcr- / I)A[j] {h)V {y + v) ■ Vf dx. 
Jr Jr Jn 

Therefore, the strong form of the boundary value problem for s is: seek s € Wi (F) 
-Kd^^s = (7 - 7<i) - / Ai [7] V (y + w) • Vfdx2 

JQ 

+ d^, A2 [7] V (j/ + w) • Vf dx2 j in F 

s(0) = s(l) = 0, 



(3.6) 



where Ai,A2 denote the representation of DA given in (2.4). We note that the 
integrals on the right hand side of (3.6) correspond to integration in X2 (vertical) 
direction. 

Finally, (13. 2b) implies 



{XuU + s, u - '«)i2(r)xL2(r) 



>0, 



Vu e Uad- 



(3.7) 



To summarize, the solution (7, y, tZ, f , s) to the first order optimality system ( |3.2[ ) 
satisfies (2.3b|, (3.51, (3.6) and (3.7). We stress that the formal approach presented 
in this section is very systematic and highly useful even though it is not clear at the 
moment how to show the existence and (local) uniqueness of the optimal control u. 
A rigorous analysis will be developed in the next two sections. 

4. The Control to State Map G^. Let Gy denote the nonlinear map 



Gy '■ Uad 

u 



{liV) 



(4.1) 



where := (F) x (fl), {j,y) solves ( p^Sbf , and the subscript on de 



notes dependence on a fixed and non-trivial v G (fl). Our goal is to show the 
existence of a control, derive the first order necessary and second order sufficient op- 
timality conditions within the realm of a rigorous mathematical framework. The first 



order optimality conditions requires to show that Gv is Frechet difFerentiable ( subsec 



tion 4.3 1 and the second order conditions require G^ to be twice Frechet difFerentiable 



(subsection 4.4 1 



The steps described above are standard for PDE-constrained optimization in fixed 
domains [1^, but our analysis for the linearized curvature OC-FBP is novel. The 
novelty resides in the highly nonlinear structure of the underlying FBP, which is 
posed in a pair of Banach spaces one being non-reflexive, and yet we deal with minimal 
regularity. A number of other control problems for FBP fall under a similar functional 
framework [T] [2], but their theory is not as complete and conclusive as ours. This 
appears to be an area of intense current research. 

The first step in this voyage is to show that there exists a unique weak solution 
to (2.3b), which implies that Gy is a well defined one-to-one nonlinear operator. In 
fact, it is known [16] that for u = and v small, a fixed point argument asserts the 
existence and uniqueness of a weak solution (7, y) in W"^ to ( 2.3b ). We further extend 
this analysis to the case where u 0. This gives us an admissible ball Uad C (F) 



where we can show the existence of solution to ( 2.3b I. Furthermore, this set guarantees 



that the state constraints in (1.2c) are always satisfied, thereby leading to a standard 



control constrained problem without state constraints. 

4.1. Well-posedness of the State System (2.3b). The weak form of the 



system (2.3b) is: find (7,?/) G W such that 



Bn [y + v,z;A [7]] =0 Vz e {ft) 

Brh,C]+Bn[y + v,EC;A[j]] ^ {u,Ow-\r).wliT) VC G M^i^r), 

where Br ■ (F) x (F) ^ M, Bn ■ {^) x (n) E are defined by 



(4.2) 



Bn [y, z; A [7]] [ A [7] Vy • Vz dx. 



(4.3) 



Furthermore, E : (F) — )■ (ft), q < 2 denotes a continuous extension such that 
E(\r = C, = (cf. [m Lemma 2]). In particular, this implies that there exists 

a constant Ge > 1 such that 



\EC\wim < CE\C\wnr) ' e (F). 
Moreover, for u E Uad C (F) we can write 

,0 



U,d:riC 

/ L==(r)xLi(r) 



(4.4) 



(4.5) 



In view of the regularity of u and C, the duality pairing in (4.5 ) reduces to uQ. This 
also enables us to deduce that for u E L"^ (F) 



l"llw'oo'(r) 



< 



<\\u\\ 



L~(r) 



L2(r) 



(4.6) 
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We will make use of these two facts repeatedly throughout the rest of the paper. 
Proposition 4.1 (inf-sup conditions). The following conditions hold for the 



bilinear forms By ■] and So [•, •; ^4 [7]] defined in (4.31 ; 
(i) Bt[-,-] is continuous and there exists a constant a > such that for every 

i&w^ (r) and s e wl (r) 



By [7,C] 



\wHT) ^ sup — , 



(4.7a) 
(4.7b) 



(ii) If ^ £ (r) satisfies (1.2c), then Bn [■r;^[7]] continuous and there exist 
constants P, Q with Q < 2 < P and /3 > 0, such that for p e (Q, P) 



I I ^ « Bn[y,z;A [7]] 



(4.8) 



Proof. For (4.7a) and (4.8) we refer to [ini Proposition 2.2-2.3] for a proof. 
For (4.7b) we proceed as follows: applying the definition of the i^-norm and the 



homogeneous Dirichlet values of s, we obtain 



\s\wl{T) 



/ |da;iS| = / Sgn(d^iS)d:riS 

Jo Jo 

^sgn(d^js) - sgn(da;,s)^ d^^s = ^Bv [C, s] , 



where C,{xi) = [^^ (^sgn(da.js) - sgn(da.js)^ € Wl^ (F). Estimate ( |4.7b[ ) follows by 

noting that ICIwi (r) — 2, and taking the sup over every C, & (F). □ 

The following lemma demonstrates how one can improve the integrability index 

of a solution to a PDE obtained by standard methods. 

Lemma 4.2 (improved integrability). Let D, be an open Lipschitz bounded domain 

of and B : (il) x Wl (17) — > M &e a continuous bilinear form. Furthermore, 

suppose that 
(i) there exists a > such that 

|xUi(a)<« sup VxGW^i(fi), (4.9) 

o^i)£Wl(n)\V\wl(^i) 

(ii) and B is continuous and coercive in W2 {^). 
Then for every F £ Wl (f2) , there exists a unique x G (il) such that 

B[xM = F{i,) for all ^ € Wl {^) and \x\wun) ^ ^WFWwim' ' (^.10) 



Proof Since (O) C Wl (17), it follows that F e (n)* C (Q)* . Lax- 
Milgram guarantees the existence and uniqueness of x G W2 (17) such that B [x, ^] = 
F{il>) for ah € (17). 
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Next, we extend [x, •] as a linear functional on (fl). To this end, let 
{i/jn} C W2 (ri) be a Cauchy sequence in the Wi (51)-norni. It immediately follows 
that {B [x, ipn]} is also Cauchy in M, i.e. 

\B [X,'0n - i'rn]] = \F{^n - i'rn)] < II ^"11 li/i (n)* I V'n " V'mliyi(n) ■ 

Finally, through density of W2 (f^) in Wi (51), not only do we obtain — t- 7/; G 
Wl (fl), but also 

S[x,V'] lim S[x,^n] - Um F(V'n) = F(V'). 



The estimate for \x\\yi (r) follows by (4.9). □ 

4.1.1. First order regularity. Now we are ready to prove that there exists a 
unique solution to (4.2) with W^-regularity. Since the system (4.2) is nonlinear we 
will obtain this result by applying the Banach fixed point theorem combined with a 
smallness assumption on a non-trivial v. To this end, we let 2 < p < P and equip the 
space = (F) x Wp {fl) with the equivalent norm 



|(7,y)||wi (l + /3C'A)|w|n/i(n)l7lwi(r) +\y\wHn) ^ 



(4.11) 



where Ca and /3 are given in (2.8) and (4.8), and define the closed (convex) ball 

:={(7,y)eWi: \y\wiin) < PCaIvI^,^^^ , < l} . (4.12) 

Furthermore, consider the operator T :My ^ defined as 

T{j,y) (Ti(7,y),r2(7,y)) = {j,y) V(7,2/) eB„, 



where 7 = Ti(7, y) £ (F) satisfies for every ( € Wi (F) 

Br [7, C] = -So [y + V, EC A [7]] + (m, C)iy^i(r)xiVii(r)' 
and y = T2{j, y) £ {n) satisfies for every z e (n) 



y + v,z;A [Ti {j,y)] 



0. 



(4.13) 



(4.14) 



(4.15) 



With these definitions at hand we proceed to find conditions under which T not only 
maps B„ into itself but is in fact a contraction in My . 

Lemma 4.3 (range of T). Let Ti and T2 he the operators defined in (4.14) and 



(|4.15|), and Ca and Ce be the constants defined in (|2.8|) and (|4.4|. Furthermore, 

(4.16) 



suppose there exists 9i e (/3Cyii/(l + /9Cyi), l) such that 

\v\w^(^n) < (1 - ^1) {o^CeCa (1 + PCa))~' 



If u £ (F) w«i/i |jM||^2(p) < Oi/oi, then the range of T is contained in My. 



Proof. Let (7, y) £ My be fixed but arbitrary. First we rely on Lemma 4.2 to show 
the well-posedness of Ti . Since it is straight-forward to check that Br is continuous 
and coercive in W2 (F), we only need to show the regularity of the forcing term in 

10 



( [4l4| ). If we define F{C) := -Bn [y + v, £(; A[-/]] + (w, C) and use ([2^, Q and 
(4.12) we find that 



F 



(C)| < Ca (\y\w^(n) +\^\wiin))\EC\w^^n) +ll"llL2(r)lClw,i(r) 

< (CeCa (1 + /3Ca)|i;U.(^) +II"IIl2(p)) ICUi(r) , (4-17) 



whence F e W^;^^ (F) and we conclude from (4.101 that 

lllwur) ^ " {<^eCa (1 + pCA)\v\w^(n) +\H\mr)) < (1 - ^^i) + - 1- 
The well-posedncss of T2 follows by Proposition |4.1| and the Banach-Necas theo- 



rem for reflexive Banach spaces ,8, Theorem 2.6]. Applying (4.8) we obtain 

\y\w^(n) < ^C'aIw liy^i(n) • 

Since (7,2/) is arbitrary, we conclude that the range of T is contained in ]B„. □ 

Definition 4.4 (admissible control set). Let 61 be as defined in Lemma 4-3 
The admissible set of controls Uad is the (nontrivial) closed ball 



Had ■■= {u e L^ (r) : ||u||i2(r) < ^i/a} • 



(4.18) 



Theorem 4.5 (T is a contraction). Let the assumptions of Lemma 4-3 hold and 
suppose further that there exists a 62 ^ (0, 1) such that 



\v\w^in) < (1 - ^2) (aCECA (1 + PCAf) 



(4.19) 



Then, the map T defined in (4.13) is a contraction in B„ with constant 1 — 92- 



Proof Consider (71, yi), (72, ^2) G B„ such that ( 71,1/1 ) ^ (72, y2)- Using ([4131 
we have that T{'ji,yi) = {^i,yi) solves (4.14| and (|4.15 l for i = 1,2. Therefore, 
combining Proposition 4.1 (i| and Lemma 4.3 with (4.191 implies 



I71 - 72lvKi,(r)< " sup Br [71 - 72, C] 



ICI 



a sup Bn[y2~yi,EC,A[ji]]+Bn[y2+v,EC,A[y2]-A[-fi]] 



CeC'a (\yi ~ y2\w^(^n) + (1 + /^C'A)l«lwpi(ji)l7i - l^lw^ir) 



= aCECA\\i7i - 72,2/1 - y2)||wi • 
Similarly, Proposition |4.1| (pl|) in conjunction with (4.15) implies 
\y\ - y2Ui(o)< sup Ba [^1 - ^2, 2; A [71]] 



(4.20) 



l^lwi(n) = l 



/3 sup So [y2 + z; ^ [72] - A [71]] 



(4.21) 



n)- 



< /3C'A|y2 + I' 1^1/^1(0)1 71 -72 1 vi/i^(r) 

< ^Ca (1 + /3Ca)|w|vvi(o)|7i -72Ui^(r) 
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Finally, ( |4.20| ) and ( |4.21[ ) imply 



(71,2/1) - (72, y2)||wi <(1 + /3C'a) klwpi(n)l7i -72lvi/^(r) 

''^livi(si) II ('^1' yi) ~ (721 2/2)1 

< (1 - ^'2) II (71, 2/1) - (72,2/2)|Li , 



where the last inequality follows from ( |4.19 1. Since 6*2 € (0, 1), T is a contraction with 
constant 1 — 02 as asserted. □ 

Corollary 4.6 (well-posedness of state system). For every u G Uad, the closed 



hall of Definition\4-4\ andv satisfying (4.16) and (4.19), there exists a unique solution 



(7, y) G W"'^ to the state equations (4.2). This further implies that Gy in (4.1) is a 



well defined, one-to-one, nonlinear operator. 

Proof. Let u S Uad be fixed but arbitrary. It now follows that T is a contraction 
in the closed convex set My (cf. Theorem 4.5) and applying the Banach fixed point 
theorem we obtain a unique (7,?/) G By such that T (7, y) = (7,2/)- In view of (4.14) 



and (4.15), this is equivalent to saying that (7,1/) is the weak solution to the FBP 
(ig, i.e. Gy{u) = (7,y). □ 



4.1.2. Enhanced Regularity of 7. 



Corollary |4.6| implies the existence and 
uniqueness of a solution (7,y) to (2.3b) with W^-regularity, provided u G Uad and v 
satisfies (4.16) and (4.19). That is, we only have one weak derivative for 7 and y. 
In the sequel we will show that the solution (7,?/) = Gy{u) is slightly more regular 
without any extra restrictions on m or u. More specifically, we will show that 



(4.22) 



The importance of this result will be evident in [subsection 5.1| where we show the 
existence of an optimal control. Despite its importance, the proof is rather simple. 
Let (7,2/) G be a weak solution to (|4.2[). Function 7 satisfies 



in the sense of distributions. If we assume, for the moment, that f e Wp (F), 
2 ^ w-i/p (p)^ rpj^ig directly implies 7 G Wp^^^^^ (F), i.e. 7 G PFp+^/' (F) 



.-1/p 



(F) as suggested. Since 



then d^^7 G VFp (F). This directly implies 7 G Wp 

as desired. Thus, it remains to show that / is in Wp 
u ^ L'^ (F), we just need to deal with the first term. 

If (/) G Wq^^ (F) then it can be seen as the trace of a function Ecf) G (fl) such 



that E(f)\r = 0, E(l>\^ = 0, and \\E(j)\\ 



With this in mind. 



(A [7] V (y + v) • I/, ^)w-'/p^r)xwy-(r) 



A [7] V (y + u) • VEc/), 



whence 



\A [7] V{y + v)- < GeCa (1 + PCaM^i 



We collect this result in the next theorem. 

The orem 4.7 (enhanced regularity). If Gy{u) — (7, y) G is the solution in 
Corollary 4.6 then 7 G VFp+^Z'' (F) n (F) without any further assumptions on u 
or V. 
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4.2. Gv is Lipschitz. The first result in showing that G„ is twice Frechet dif- 
ferentiable is to demonstrate that it is Lipschitz. 

In the interest of saving some space we will rewrite the variational system ( |4.2[ ) 
in the following form: find (7, y) e such that for every {(, z) € Wl (F) x (f2) 

Br[i.(:]+Bn[y + v,z + E(:-A[-f]] = {u,Ow-\r)>cwl(rr (4-23) 
With this new notation in place we are ready to study the Lipschitz continuity of Gy. 

Theorem 4.8 (Lipschitz continuity of Gy). If v fulfills the conditions of Corol- 
lary \4-.€\ then Gy satisfies 

\^Gy{ui) - Gy{u2)\\^^ < Lg\\ui- U2\\^2(Y) "iui,u2 eUacu (4-24) 

with constant Lq — f^{^ PGAf'\v\Y/i[Q,)- 

Proof Given mi, U2 € Uad, set (71, yi) - (72, 2/2) = Gy{ui)- Gy{u2). Using ( |4.23[ ), 
we have for every {(^, z) e Wl (F) x (17) 

Bt [71 - 72, C] + Bn [yi +v,z + EC,A [71]] 

-Bn [2/2 + t', z + £^C; ^ [72]] = ("1 - "2,C)L2(r)xL2(r)- 

Subtracting Bq, [y2 + v,z + EQ; A [71]] from both sides and rearranging terms yields 

Bt [71 - 72, C] + So [vi -y2,z + EC, A [71]] = 

+ Bn [y2 +v,z + EC;A['j2] - ^ [71]] + (^i - "2, C>L2(r)xL2(r)- 



The inf-sup estimates from Proposition 4.1 together with (ji^yi) G B„ for i = 1,2, 
imply for C = 

1 2/1 - 2/21^/^1(0) < /3C'a|2/2 +t'lw^i(a)|7i -72lwi(r) 

< 13Ga (1 + /3C^)klw^i(o)l7i - 72Ui(r) , (4.25) 

and for z = 

I71 -72lM/i(r) < ^CeCa (\yi - 2/2lw;(n) +I2/2 + w|vy^i(n)l7i -^^Iw^^r)) 
+ a\\ui - U2||^2(r) 
< uCeCa (1 + /3CA)^|f |,y^i(o)l7i - 72lvi/^(r) + "Iki - "2|li2(r) • 



Finally, in view of (4.19), we infer that 



I71 -72lM/^(r) < ^ll^ii - ^i2|lL2(r) • (4.26) 



The asserted estimate follows immediately from the definition of ||-||^i in (4.11). □ 

4.3. Gy is Prechet Differentiable. The next step towards showing the twice 
Frechet differentiability of Gy entails analyzing the well-posedness of the linear vari- 
ational system: find (7, y) G W"^ such that for every z) G (F) x (fi) 

Br [7, C] + % [(7, 2/) , z + EC; 7, v] = En {z + EQ) + (C) (4.27) 
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where 



[(7, y) , •; 7, y] Bn [2/, •; A [7]] +Bn[y + v, •; DA [7] (7)] , 

(7, y) = Gy{u) € B„ for a fixed u inUad, DA [7] (7) is given in Q, and Fn G VF^^ (rj)* 
and Fp € W{ (F) are fixed but arbitrary. 



4.3.1. Preliminary Estimates. Given tliat tfie coupled system (4.27) is linear, 
one would be inclined to use the standard Banach-Necas theorem to prove its well- 
posedness directly. We deviate from this approach and resort to the machinery already 
put in place. 

Consider the operator T : given by 

r(7,y):=(ri(7,y),r2(7,2/)) =(7,y) V(7,y)eW\ (4.28) 
where 7 = Ti (7, y) e (F) satisfies for every G (F) 

Br [7, C] = [(7, y) , ^^C; A [7]] + Fn {EQ + Fr (C) , (4.29) 
and y (7, y) G Wp (il) satisfies for every z G (fi) 

Bn [%z;A [7]] = -Bn [y + v,z;BA[j] (Ti (7, y))] + Fn (z) . (4.30) 



We point out that any fixed point of T is a also a solution to (4.27). To infer the 
existence of a fixed point we exploit the linear structure of (4.27). Therefore, it suffices 
to show the well-poscdness of the intermediate operators Ti and T2, and to show that 
T is a contraction . 

Lemma 4.9 (well-posedness of Ti and T2). Let Ti, T2 he the operators defined in 
(4.29) and (4.30) with (7, y) G The following holds 



(i) for every (7,?/) G W""^, there exists a unique 7 = Ti (7,?/) satisfying (4.29) and 
\l\wl^{T) ^ " {CECA\\{l,y)\\^i + C£;||Fj2||^yi(s-j). + ll-^r || h^^i (p)' ) , 



(ii) for every (7,2/) G W"'^, there exists a unique y ~ T2(7,i/) satisfying (4.30) and 
\y\w}{n) < /^C'a (1 + ^CA)|w|vy^i(n)l7lvi/4(r) + l^WP^iWw^Hny ■ 



Proof. To prove ([i|) we proceed as in Lemma 4.3 It suffices to check that the 
right hand side RHS(C) of ( |4.29[ ) is in Wl (F)*, namely 

|RHS(C)| < (cECA{\y\y^^^a) + i^ + PCA)\v\^^^n)h\wi^i^)) 



CE\\Fn\ 



'+\\Fr\ 



4.2 



with the coercivity of Br in W2 (F), 



The desired estimate follows from Lemma 

and the definition of |j-|j^i in (4.11). 

Estimate (|ii| is a straightforward application of the Banach-Necas theorem [8] . □ 
Theorem 4.10 (T is a contraction). Let ( |4.19[ ) hold for some 62 G (0,1). The 

operator T defined in ( 4.28[ ) is a contraction inW^ with constant 1 — 02- 
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Proof. We proceed in a similar fashion to Theorem |4.5| Consider not identical 
(71, yi) and (72,2/2) in W\ and use ( |4.28| ) to write (ji^yi) = T{-fi,yi) for i = 1,2. 
Applying Lemma 4.9 (i|, we obtain 

I71 -72lw^(r) < aCijCA II (71 - 72,2/1 -2/2)||wi • 
Similarly, Lemma |4.9| ^ implies 

|yi - hlw^n) < (1 + /3C^)kliyi(o)l7i " 72lvi/^(r) ■ 



Lastly, the upper bound (4.19) on v yields 

||(7i -72,yi - y2)||wi < (1 + /3C'A)^|t'lvy^i(n)|7i -72lvKi(r) 
< (1 - ^2) II (71 - 72,2/1 - 2/2)11^1 • 
Hence, T is a contraction with constant 1 — 6*2, as asserted. □ 



Corollary 4.11 ( well-posedness of the linear system ( 4.27 ) ) . Under the assump- 
tions of Theorem 4-10, there exists a unique solution (7,2/) € to the variational 



equation (4.271 and the following estimates hold 
a 



< 



II < ^ 
\y\wi(n) ^ 

Therefore 

||(7,2/)|Li ^ 



Ce (1 + /3C'^)||^o|lvF,i(f2)' +ll-^r|lvVi(r)" 
(\\Fn\\wi{nr + ^Ga{1 + PCA)\v\^i^n)\\Fr\\wliry 



(4.31) 
(4.32) 



^ [aCE (1 + /3CA)>lH.;(n) + /?) WMw^iny 



Proof. Existence and uniqueness follows from Theorem |4.10[ As far as the esti- 
mates, we will only derive (4.31) since the other two are merely a consequence of it. 
To this end we apply Lemma 4.9 and the upper bound (4.19) for v to get 

blw^iD < o^CeCa ((1 + PCa)\ ^Iwi(n)l7lw^(r) j 
+ " {^^EWFnWwKny +ll^r|lH^i(r)*) 

< aCECA{l + f3CAf\v\^ri{n)\l\wwc) 

+ aCE{l + pCA)\\Fn\\w^^^nr+a\\FT\\wl(T)' ■ 

< (1 - ^2)\l\w}^^T) + ^Ce (1 + PCAWnWwiiny + "ll^r|la.i(r)* • 



The estimate (4.31) follows immediately. □ 

4.3.2. The First Order Prechet Derivative. In this section we will prove the 
first order differentiability of the control to state map G^. Since the set Uad defined 
in (4.18) is not open we need to define a proper set of admissible directions. 

Definition 4.12. (admissible directions) Given u € Uad, the convex set C (u) 
comprises of all directions h G (F) such that u + h G Uad, i-G., 



C{u) := i^hGL^{T) : u + hGUad] 
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Theorem 4.13 (the Frechet derivative of Gy). The control to state map Gy : 
i^ad — >■ W"'^ is Frechet differentiable. Its first order derivative (7, y) := Gy (u) h in 

at an arbitrary u in Uad in the direction h in C (u) satisfies the linear variational 
system (4.27) with = and Fr (C) — /p ^C^ namely 

Br [7, C] + T^n [(7, y),z + EC; l^v]^ V (C, z) e Wl [T) x (n). (4.33) 



Moreover, the following estimate holds 



|(7,y)| 



(4.34) 



Proof. The derivation of (4.33) is tedious but straightforward, so we skip it. The 



estimate (4.34) foUows from CoroUary 4.11 



We turn our focus to proving that G'y is the Frechet derivative of Gy. To this 
end, we must show that the remainder operator CHg^ , defined as 

y{a^[u,h]:=Gy{u + h)-Gy{u)~G'y{u)h yu e Uad, h e C {u), (4.35) 

satisfies for all u £ Uad 

||^G„ [u,h]\\ 



lim 

?ieC(M):||/l||^2, 



\h\\ 



0. 



Since we do not have direct access to \\^g^. strategy of the proof is to 

first show that 9^g„ [m, h] satisfies ( |4.27| for some Fn € {n)* and Fr = 0. Next, 
owing to the estimates in Corollary 4.11[ it suffices to check that 



lim 



0. 



h<£C(uy.\\h\\^2fr)^a ||/i||i2(r) 

To avoid any ambiguity we adopt the following notation in this proof, 

(7(?2), y{u)) := Gy (u) , {^{u + h),y{u + h)) := Gy{u + h) 

{lu{u)h, yu{u)h) := Gy{u)h, {S-y, Sy) = (<K^ [u, h] , ^Ry [u, h]) 9^g„ (m, h) , 

whence 

^7 = 7 (t2 + /i) - 7 (u) - 7„ {u)h 5y^y[u + h)-y {u) - ?/„ [u) h. 



According to the remainder definition in (4.35) we start by combining (4.2) for 
Gy {u + h) and Gy {u) with to obtain for every (C, z) in Wl (F) x Wl (fl) 

= Br [jiu + h) - j{u) - 7„(u)/i, C] - % {lu{u)h,yu{u)h) , z + EC,\-f{u),y{u) 



Br. 



y{u + h) +v,z + E(;A [7(1* + h)] 



Be 



y{u) + v,z + EC;A[j{u)] 



Adding and subtracting Vn (^^{u + h) ~ -f{u),y{u + h) — y{u)) , z + EQ; j{u),y{u) 
to the previous equation and utilizing the definition of 6^ and 5y above, yields for 
every (C, z) in Wl (F) x (n) 

Bn [<57, C] + % [{Si, 5y),z + EC; liu), y{u)] = Fn{z + EC), 
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where 



Fni-) = Bn y{u + h) + v,-;A [-fiu)] - A [7(7/ + h)] 



Br 



y{u) + V, •; [j{u)] {j{u + h) ~ 7(u)) 



The fact that Fq is in {fl)* follows from the continuity of [w, •; V] with|w|^i 



and II y I 



L°°(n) 



bounded uniformly (c.f. (2.8)). Our last step is to add and subtract 



y{u + h)+v, •; DA [7(72)] (7(72 + h) - j{u)} 



to F^, employ the definition of the 



Bn 

remainder (|2.5|) and the Lipschitz estimates (|4.25|) and (|4.26|) to obtain 



lim 

heC(«):it/il|^2(r)^0 



d\A [7(u), 7(w + h) - 7(w)] 



L°=(0) 



as well as 



\^{u + h) — j{u) 



+ CA\y{u + h) - y{u)\ 

= o(ll'i|lL2(r)^ 
This concludes the proof. □ 

4.4. The Second Order Prechet Derivative. The final result of this section 
is to show that G„(u) is twice Frechet differentiable with respect to u. The technique 
is similar to the one in subsections |4.2| and |4.3| In fact, proceeding as in Theorem |4.8| 
we get the following. 

Proposition 4.14 (Lipschitz continuity of G^). There exists a constant Lc > 0, 
such that for every ui , 1*2 G Uad 

\\G'^{ui)h~G'Ju2)h\ 



sup 

0^hi£C{ui)nC{u2) 



< Lg'\\ui - U2I 



lL2(r) 



LHT) ■ 



(4.36) 



Theorem 4.15 (the Frechet derivative of {u)h). The control to state map 
Gv '■ l^ad ^ W""^ is twice Frechet differentiable. Its second order derivative (7, y) := 
G" {u) hih2 in at an arbitrary u G Uad in the direction (hi, /12) G C (w) x C (u) 
satisfies the linear variational system (4.27), namely for every {C,z) in (T) x 

Br [7, C] + % [(7, y) , ^ + EC; 7, y] =Fn{z + EQ , (4.37) 
with Ffi £ (^)* given by 

Fni-) :--6o[z/i,-;DA[7] (72)] 

r 9 1 (4-38) 

- So [2/2, •; Dv4 [7] (71)] - + t-, •; EPA [7] (71, 72)J , 

and {"1i,yi) '■— G^ {u) hi, for i = 1,2. Moreover, the following estimates hold 

l7U^(r) < ^CeCa (1 + 2PCa) (1 + /3G^)'bl wi(n)ll^illL2(r)ll^2|lL2(r) ' (4-39) 
a" 



\y\w^(n) < -fpPCA (1 + 2/3Ga) (1 + PCa)\v\ 
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(4.40) 



4.11 



with 



Proof. We skip the derivation of (4.37) because it is tedious but straightforward. 
The estimates for |7|^yi (-pj and |y|i^i(f2) are a consequence of Corohary 
Fr = after estimating (4.38), namely 

\\Pn\\w^(_n)- < Ca (|yilwpi(n)l72lvy^(r) +\y2\w}in)hi\w^{r)) 

l7i| 



< 2 



- pCiil + f3CA)\v\y^.^n)\\hi\\ L2(r)ll'*2|lL2(r) 



CAil + fiCA)\v\wHn)\\hi\\ L2(r)ll^2|li2(r) 
Ca (1 + 213Ca) (1 + l3CA)\v\^.^n)\\hi\\LHr)\\h2\\mr) ■ 



where we have used (4.31 1-(4.32 1 with Fn = for {jiTyi) along with (4.6 1. 

The strategy for showing second order Frechet differentiability of Gy is the same 
as in [Theorem 4.131 wc first show that the remainder 



(^7, Sy) := {u + h^) h, - G'^ {u) - G'l {u) h^h^ 



(4.41) 



satisfies the linear variational system in (4.27) for a suitable right-hand side 6Fn G 
(i^)*, and prove that 



sup 

o=^hiec{u) 



= oi\\h2\ 



i^illL2(r) 

As a tradeoff between clarity and space we denote uq ~ u, Ui — u + hi, and 

(Tu K) hj,yu {ui) hj) (ui) hj 

{luu (uo) hih2,yuu (wo) /11/12) G" (uq) /11/12, 



(4.42) 



for i,i — 1, 2, whence 

h = In ("2) hi - 7„ (wo) hi - 7„„ (uo) hih2, 
Sy = y„ (U2) hi - yu (uq) hi - (wq) /ii^2- 

According to the remainder definition in (4.41) we start by combining ( 4.33[ ) for 
G'„ (ua) hi and G'„ (uq) /ii with ( |437| to obtain for every ((, ^) in M^ii (F) x Wf{n) 

-Fn {z + EC) = Br [Sj, C] + % [G'^iu2)hi,z + EC; G^iu^)] 

- Vn [G'^{uo)hi + Gl{uo)hih2, z + EC; G„(uo)] , 



where — -F'n(-) = defined in (4.38). Further manipulation, based on 

adding to both sides the following two additional terms, 

F^{z + EC)^Va {G\Xu2)hi,z^EC;Gy{uQ)\ , 
F5 (2 + EC) - -Vn {G',{u2)hi,z + £;C; G,(u2)] 
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leads to 



T-2 - Br 



T3=Bn 



n - Br 



Te=Br 



Br [Sj, C] + ■Dn [(^7, Sy),z + EC, 7(^0), y("o)] = SFn {z + EC) , 

where SFq ~ Si=i -^i clearly in {^)* ■ To create additional cancellations we 
further decompose SFq — Y^^^^i'^i follows: 

Ti = Bn yu{uo)hi, ■■,'DA [7(^0)] (7('«2) - 7(^0) - 7«(wo)/i2) 
yu{u2)hi - yu{uQ)hi, •; DA [7(^0)] (7('"2) - 7("o)) 
yu{u2)hu •; A [7(^2)] - A [7(^0)] - DA [7(1*0)] (7(^2) - 7(^0)) 
yu{uo)h2, •; DA [7(^0)] (7„(u2)/ii - lu{uQ)h 
y{u2) - y{un) - Vu{u(3)h2, •; DA [7(^0)] {lu{u2)h 

DA [7(^0)]) (7 

-D?A [7(1/0)] (7«(«2)^i,7("2) - 7(wo)) 

?7 So [^(""2) + •; D?A [7(1*0)] (7«("2)/ii, 7(^*2) - 7("o) - ^uiuo)h 

y{u2) + V, sD^A [7(^0)] (7„(u2)/ii - 7M(^o)/ii, 7«(wo)/j2 

J/(w2) - y{uo), -l^A [7(1*0)] (7«("o)/ii,7«("o)^2 
We now estimate each of these terms separately and show 



Ts - Bn 
Tg - Bn 



II II 

0#/iiGC(«o) ll"-lllL2(r) 



lL2(r) 



(4.43) 



which obviously imply (4.42 1. 
• Term Ti : Since 



\\Ti\\wi^ny < CaIvu (uo) hi\^,^^^\-f {u2) - j (uo) ~ (wo)/i2|^j^(r) , 



the estimate (4.34), together with 

I7 ("2) - 7 (wo) - 7m (wo) /J2 



vi/^(r) 



= O (l|/l2|lL2(r) 



implies ( |4.43 |. 
Term T2: Since 



||72|lH'i(n)- < C'AlyM (U2) /ii - y„ (uo) /''i|vi/i(n)b ("2) - 7 K)lvKi^(r) 



it suffices to recall the Lipschitz properties (4.24) of Gy, and (4.36) of to deduce 
( 4A3| . 

Term T3: Invoking the Frechet differentiability (2.5) of A, and the Lipschitz prop- 
erty (4.261 of 7(uo) we infer that 

A[7(u2)] -A[7(ko)] -DA[7(iio)] (7("2)-7("o)) ^ 

(Si) 

= o(|7(u2)-7("o)|^i^(r)) = o(ll'*2|lL2(r)) • 
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This, in conjuction with (^2) 
• Term T4: In view of the Lipschitz property (4.36) of G'^ 



yields (4.43) 



|7«(u2)/ii -7«K)/ii|vi/i^(r) ^ll^illL2(r)ll^2|lL2(r) 



property ( |4.43[ ) follows from \yu (uq) ^2|pyi(j^) ^ ll^2|lL2(r)- 



• Term T5: Since y{u) is Frechet differentiable according to Theorem 4.13 namely 

^l|/»illL2(n- 



the bomid (4.43) is a consequence of |7„ (1*2) j-p-) 
• Term Tq: We recall the second order Frechet differentiability (|2.7| oi A , and the 



Lipschitz continuity (4.24 1 of G^, to write 



[7 (m2)] {lu{u2)hi) - DA [7 (Uo)] (7n(w2)/ll) 



■D^-^ [7(^0)] (7«("2)^i,7("2) -l{uo)) 



= \\hi\ 



Since |y(M2) 



L2(r) 
<l 



' (I7 (W2) - 7 (wo) 



L~(r) 



= \\hi\ 



L2(r) ' 



lL2(r) 



, this implies (4.43 ) 



Term T7: We proceed as with Tg, now appealing to (2.8) and the Frechet differen- 



tiability of 7 at uq (Theorem 4.13), to obtain 
|7«(u2)/ii|i^^(P)|7(w2) - 7(wo) - lu{uo)h2 



= \\hi\ 



0(11 



whence (4.43). 



• Term Tg: We employ the Lipschitz property (4.36) of in to write 

\lu{u2)hi - 7«(uo)/Ji|vi/^(r) ^ll'iillL2(r)||/i2|li2(r) • 

<l 



The desired bound (4.43) follows from 7„ (uq) /12 



lL2(r)- 



Term Tg: We use the Lipschitz property (4.24) of G^ 
\y{u2) - y{uo) 



together with (mo) (r) ^ Inill r2fr\ to deduce (|4.43|) 



•'lw'i,(r) 

Altogether, this concludes the proof. □ 



<\\h9. 



L2(r) 



5. Optimal Control. Let us summarize what we have accomplished so far. We 



have formally derived the first order necessary optimality conditions in section 3 If 
Gv denotes the control to state map, we have proved in [section 4| that G^ is well 
posed, i.e., there exists a unique weak solution to the state equations (2.3b) for every 



u e Uad in (4.18), and v satisfying (4.16) and (4.19). As a crucial step forward we 



have shown that Gy is twice Frechet differentiable. 

This background work puts us in the position to show the existence and uniqueness 



of the optimal control u solving the OC-FBP in (2.3a)-(2.3b). We will achieve this 



result in three stages. We first show the existence of u in Theorem |5.1| of [subsection 5.1[ 
We next derive the first order necessary optimality conditions and show the existence 
of solution to the adjoint equations in [subsection 5.2[ Finally in [subsection 5.3[ we 
end this voyage by proving the second order sufhcient conditions for the control u. 
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5.1. Existence of Optimal Control. In order to show the existence of a 
solution to our optimal control problem we must first rewrite the cost functional 
J : X Uad from (2.3a) in its reduced form. This is accomplished by utilizing 
the control to state map G„ from Section [4] as follows: 

Jiu) JiG4u),u) = J{j,y,u) = ^i(G„(u)) + J2(m), 

with 

1„ ..9 X 



;ll7-7d|lL2(r) + -^Wv + ^ - Vdh^n) > ^2(w) 



2 

L2(r) 



Now an equivalent minimization problem to (2.3a) is 



min J^{u) 



(5.1) 



Theorem 5.1. For every v satisfying (4.16 
control u G Uad minimizing the cost functional 



1+1/9 



(4.191, there exists an optimal 
2.3a I with optimal state (7,?/) € 



(r) n (r)j x (ri) which solves the free boundary problem (2.3b) and 

satisfies the state constraint (1.2c). 

Proof. In order to show the existence of an optimal control we use the direct 
method of the calculus of variations. We first note that the cost functional J' in (5.1 1 
is bounded below by zero, whence j — iniueUad J^i'^) finite. We thus construct a 
minimizing sequence {un}„(=N such that 



J 



lini J{ur, 



By Definition 4.4 Uad is nonempty, closed, bounded and convex in (F), thus weakly 
sequentially compact. Consequently, we can extract a weakly convergent subsequence 
{wrifclfeeN C (r), i.e. 



inL^(r), ueUad- 



Here u is our optimal control candidate. 

Hencefor th, wc drop the subindex k when extracting subsequences. According 

to Corollary 4.6 and (|4^, we let Gy{u„) = (7„,2/„) G (Wp^^^'' (T) n W^^ (F)) x 
Wp (fl) denote the unique state corresponding to w„, thereby solving the free boundary 
problem ( |2.3b ) and satisfying the state constraint ( 1.2c). Since W^^^^'' {r)C\W^ (F) is 
compactly embedded into (F) , using Rellich-Kondrachov theorem we can extract 
a strongly convergent subsequence {7,i}„eN C W}^ (F), i.e. 



in Wl (F), and 



in Wl (n). 



Note that the limit pair (7, y) is the state corresponding to the control u. This results 
from replacing (7,2/) with (7„,2/n) in the variational equation (4.2) taking the limit, 
and making use of the embedding i^(F) C (F). 

Finally, using the fact that J2 (u) is continuous in and convex, together 
with the the strong convergence {'jn,yn) {j,y) in (F) x (il) due to Rellich- 
Kondrachov theorem, it follows that J7 is weakly lower semicontinuous, whence 

inf ^(u) =liminf (^i(G„K)) + J2K)) > Ji{G,iu)) + J2{u) ^ J (u) . 

u^Uad n— >-oo 



This concludes the proof. □ 
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5.2. First Order Necessary Condition. We start with a classical result [19]. 

Lemma 5.2 (variational inequality). Ifu£ Uad denotes an optimal control, given 
by Theorem \5.1\ then the first order necessary optimality condition satisfied by u is 



(J'iu), 



(5.2) 



We will show that the variational inequality (5.2) is the same as (3.7) and in the 
process we will prove that (3.51 and (3.6) are the correct adjoint equations, thereby 
furnishing a rigorous derivation of the formal results in Section [Sj 

Theorem 5.3 (first order conditions). If u (z Uad denotes an optimal control of 
OC-FBP, then the first order necessary optimality conditions are given by (3.5), (3.6| 
and (3.7). 

Proof. We can infer that J is Frechet differentiable by recalling from Theo- 
rems |4?T3] and |4?T5] that Gt, is twice differentiable and that Ji is quadratic. In fact, 
the Frechet derivative of in (5.1 1 at w in a direction h is 

J'{u)h = j[{G.,{u))G'vi^)h + J2{u)h = j[{G.,{u)){lu,Vu) + J2{u)h, 



(5.3) 



where {'fu, yu) = Gy{u)h satisfies (4.331 and 

J'{u)h ^ Xy{y + v^ yd, yu)L2(fj)xi2(o) 

+ (7-7d,7u)L2(r)xL2(r) + ^u{u,h) ^2^^)xL^^)■ 
lntroducmg the adjoint states (f, s) e (51) x (F), which satisfy the system 
(3.5)-(3.6l in weak form, and noting that h e i^(F), we obtain 

J'{u)h = Br hu,s] +Vn [hu,yu) ,r;7,y] + A„(u, /i)i2(r)xL2(r)- 



Utilizing (4.33) with C — s and z — r, we arrive at 

J'{u)h = (s + A„?2,/i)^2(r)xL2(r) [{lu,yu) ,r - Es;j,y] 



Since the Dirichlet condition r|r = s implies r — Es G (O), (4.33) with C = and 
z e Wl {Vl) yields [(7„, ?/„) , z; 7, y] = 0, whence 



J'{u)h = (s + A„w, h)j^2 



L2(r)xL2(r)- 



In view of (5.2 1, this coincides with (3.7) for h = u — u admissible. □ 

5.2.1. Well-posedness of the Adjoint System. Before we dwell upon the 
second order sufficient optimality conditions we put together the last piece of the 
puzzle: the well-posedness of the adjoint system (3.5) and (3.6). This will be done 
using a contraction argument in Banach spaces, assuming that we have a solution 



(7, y) € My to the state equations in (2.3b) satisfying Proposition 4.1 and Lemma 4.3 
Let V := |r e {il) : r\r G (F), r|s = o|, and the operator Ti : V H 
Wi (F) be defined as s = Ti (r) where s satisfies for every ( £ (F) 
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w'iCrj.Wf ^(r) 



(5.4) 



with 



and 



/ Ai [7]Vy- Vrdxa, 
Jo 



/i[-;7,y,''] 



/ ^2[7]Vy- Vrdx2. 



(5.5) 



Given s € (F), let : W^^i (F) ^ V be the operator defined as¥ ^ T2(s) = i + Es 
with ^ e W"^^ (ft) satisfying 



-4 [7]] = Xy{z,y + v - yd) 



(0)xL9(n) 



Vz e (f]). (5.6) 



Lemma 5.4 (ranges of Ti and T2). Let Ti, T2 &e defined in (5.4) (5.6). // 
(7,y) S B^,, defined in (4.12), then for every 
(i) r eY, the solution s = T'i(r) to (|5.4[) satisfies 



\s\wiir) < " (II7 - 7d|lLi(r) + 2/3C'lklvi/pi(o)l'"U,i(f^) 



(^zzj s G (-T); the solution r = T2 (s) to (5.6) satisfies 



\r\w^{n) < (\\y + - ydhiiu) + ^'eCa\s\wi 



(r) 



Proof. Using (4.8) of Proposition 4.1 and applying Banach-Necas theorem [8], 
there exists a unique solution rto (5.6). Estimate ([ii| follows from (4.8) and (2.8), as 
well as the Poincare inequality II^^H^pfQ-j < klvyi(f2) the unit square. 

In order to show the existence of solution to ( |5.4[ ) we note that we are looking for 
an absolutely continuous function on P = (0, 1) with zero Dirichlet values. Therefore, 
by the characterization of such functions in M p[5J Theorem 5.14], there exists a 
5 e (0, 1) such that 

s{x,) = ^ ' (g(i) - ^ .g(r) dr) dt Vx, € (0, 1) 

because si(0) = si(l) = 0. The variational equation (5.4) is satisfied by swith 



git) 



1 



fi[t-l,y,r]+ fo[T;j,y,r]di 



and /o,/i defined in (5.5). 

It remains to check that /J fodr and /i arc in (0, 1). This follows by applying 
Fubini's theorem and Holder's inequality. We consider first 

/ |/o[a;i;-]|dxi <||7-7d||^i((,i) + / |Ai [7] Vy • Vr| dx 
Jo Jn 

< II7 - 7<i|lLi(04) + CA\y\w^(n)Mw^{n) 
<\h - Idh^oA) + ^'^A\v\w}(n)Mw^{n) - 
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because < l3CA\v\^ri(^Qy Similarly, we obtain an estimate for /i 



This implies Finally, the uniqueness of s follows from the estimate in (4.7b). □ 

Theorem 5.5 (existence of (3.5) and (3.6)). Under the assumptions of 
Lemma 5.4' the operator T ~ T2 o Ti : V — > V is a contraction with constant 
l3CA/{l + f3CA). 

Proof. The proof is similar to the one in Theorem |4.5[ therefore we will be brief. 
Consider ri, r2 G V such that ri 7^ r2 and let 's.^ = Ti{ri),ri = T2(si), where Si,rj 
solve (5.4) and (5.6) for i = 1,2. Then Proposition |4.1| (p| and Lemma 5.4 (!]) imply 

|si - S2lvvi(r) < a sup Br [C, si - S2] 

Klw^(r) = l 

= a sup (C,f [7,y,ri] - f [7,y,r-2]) 



ICIv 



In addition, since 



\r1-r2\wiin) 



£1 + Ell -h- Es2 



< C 



and Proposition 4.1 (ii| implies 



< /3 sup Bn 

9^ ' l^lw^in) — ^ 



e\Si - S2lvyi(r) 



z,h-£2;A[j] 



= 13 sup Bn[z,Es2- Esi;A[y]] , 

we deduce that 

\ri - r2\w^^n) < C'b (1 + PCa)\si - S2lvi/i(r) 

< aCECA (1 + I3Ca) /3C'A|w|i4.i(n)ki - ^2Ui(n) 



Invoking (4.16) along with 9i > f3CA/{^ + PCa) we obtain 

\ri-r2\wi(n) < (1 - ^'i) /^C^lri - r2|^_i(f^) < ^ _^ ki - ^^2 1 vy^i (o) 



Therefore T = T2oTi 



is a contraction in V. □ 



5.3. Second order sufflcient condition. The final step is to prove the second 
order sufficient condition for the optimal control u found earlier, which in turn guar- 
antees that u is locally unique. This hinges on an additional condition on Iwl^^i/i^fj) 
depending on the parameter 6*3 given by 



= ^ {Ca (1 + 13Ca) (1 + 2(3Ca)) ' 

• Lce (1 + PCa) (1 +||7<i|lL2(r)) + +\\yd\\ 



L2(0) 



(5.7) 
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Theorem 5.6 (second order sufficient conditions). IfOi, 62 satisfy (4.16), (4.19) 
and in addition 



< 



then 



(5.8) 



(5.9) 



Proof. Since Gy is twice Frechet differentiable, according to Theorem |4.15[ we 
can write the second order Frechet derivative of J from (5.1 ) at m in the direction 

as 

J"{u)h^ = Jl'{Gy{u)){G',iu)hf + Jl{G.,{u))G'^{u)h^ + Ji'{u)h''. (5.10) 



By recaUing that {■ju,yu) — G'y{u)h solves (4.33) and {■^umUuu) = G'l{u)h? solves 



(4.37), we can write 



where 



J"{u)h' = Ji"(G,(S))(7„,y„)2 + J[{Gy{u)){^uu,yuu) + J2{u)h\ 

Ji{Gv{u)){-iu,Vuf =ll7«lli2(r) +Aj,||y„||^2(o), 

Ji'(Gt,(u))(7„„,y„„) = / (7-7d)7«M + Ay / {y + v - yd)yuu, 
Jr Jn 

This yields 

J"{u)h^ > A„||/i||^2(r) +11 7tt|lL2(r) + Ay 1 1 j/u 1 1^2(^2) 

-II7- 7<i|lL2(r)l7««Ui,(r) -KWy + 'v- yd\\L^(^n)\yuu\w}(n) , 

because of the Poincare inequalities ||2:||i2(f2) — klwi(O) — ^"^^ ^^"^ square 

n and||C||i2(r) <IClM/i(r) <\C\w}{r) t^e unit interval F. 
Estimates ( |4.39| ) and ( |4.40[ ) in Theorem [4. 15| imply 



J"{u)h' > A„||/i||i2(p) - (1 + (3Ca) (1 + 2|3CA)\v\^^^^^ 

■ (aGE (1 + I3Ga)\\i - IdW L2(r) +/3Aj;||y + v- y^W ^2(0)) I|/i|li2(r) 
Furthermore, from (7, y) G By we obtain 

ll7-7d|lL2(r) <||7llL2(r) +ll7d|lL2(r) < 1 +\hd\\L2(n) > 

and 

\\y + v-yd\\L2(^^j <|y + v|vi/^i(o) +Ily<i|li2(s^) 
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with 



\y + v\y^ifm < {I + ^3CA)\v\y^,^ 



< 



1 - I 



In view of the definition of ^3 above, this in turn impUes 



|2 



Therefore, the smallness condition (5.8 1 on 1; yields (5.9 1. □ 

Corollary 5.7 (quadratic growth). Under the assumptions of Theorem 5.6 
there exist 6 > such that for all h £ C (u) wit/i ||/i||^2(p) < we have 



J{u + h)>j(u) + ^\\h\\l.^^y 



(5.11) 



Proof We refer to fW, Theorem 4.23]. □ 

Corohary |5 . 7| imphes that there exists a unique local minimum u solution to our 



— u||^2(p) yuEu + C{u). (5.12) 



OC-FBP. Moreover, (5.11) is equivalent to 



{J'{u)^ J'{u),u~u) 



L2(r)xL2(r) 



T 
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